close all
clear all

%% Switch on/off (true/false) of special features
p.dostructboost = false;
p.premrsgfcy_is_rw = false;

setparam;

%% Model
m = model('model.model','assign',p);

%% Check steady state
[flag,discrep,eqtn] = chksstate(m,'error',false,'tolerance',1e-8);

if ~flag
  error('Equation fails to hold in steady state: "%s"\n', eqtn{:});
end

%% Model solution
m = solve(m);